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ABSTRACT. The problem is: Let Hn = 7 be the n-th harmonic number. Show, for each 
n > 1, that 

^d<Hn + exp{Hn) log{Hn), 

d\n 

with equahty only for n = 1. 
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1. Introduction 

We consider the following problem. 

n 

Problem E. Let Hn = X] ^- Show that, for each n> 1, 

J2d<Hn + eMHn)^ogiHr,), (1.1) 

d\n 

with equality only for n = 1. 

The function cr{n) = X]^|„ d is the sum of divisors function, so for example cr(6) = 12. The 
number Hn is called the n-th harmonic number by Knuth, Graham and Patashnik [^, Sect. 6.3], 
who detail various properties of harmonic numbers. 

The in Problem E might stand for either 'easy' or 'elementary'. Perhaps 'iJ' for 'hard' 
would be a better letter to use, since our object is to show the following equivalence. 

Theorem 1.1 Problem E is equivalent to the Riemann hypothesis. 



The Riemann hypothesis, stated by Riemann in 1859, concerns the complex zeros of the 
Riemann zeta function. The Riemann zeta function C,{s) is defined by the Dirichlet series 



as) = E 



n 

n=l 



which converges for 3f?(s) > 1, and it has an analytic continuation to the complex plane with one 
singularity, a simple pole with residue 1 at s = 1. The Riemann hypothesis states that the nonreal 
zeros of the Riemann zeta function C,{s) all lie on the line 3?(,s) = |. One reason for the great 
interest in the Riemann hypothesis, regarded by many as the most important unsolved problem 



in pure mathematics, is its connection with the distribution of prime numbers, described below. 
More significantly, the Riemann hypothesis is a special case of questions concerning generalizations 
of the zeta function (L-functions) and their connections with problems in number theory, algebraic 
geometry, topology, representation theory and perhaps even physics, see Berry and Keating Q, 
Katz and Sarnak [| and Murty[|ll|. 

The connection of the Riemann hypothesis with prime numbers was the original question studied 
by Riemann Let 7r(x) count the number of primes p with 1 < p < x. C. F. Gauss noted 

empirically that 7r(x) is well approximated by the logarithmic integral 



Li{x) 

which itself satisfies 



dt 
logt' 



Li{x) = -^+0{- ^ 



logx (logx)^ 

The Riemann hypothesis is equivalent to the assertion that for each e > there is a positive 
constant such that 

\tt{x) - Li{x)\ < Qx^/2+. 

for all X > 2, see Edwards p. 90]. The force of the Riemann hypothesis lies in the small size of 
the error term. The strongest form known of the Prime Number Theorem with error term asserts 
that 

\ir{x) - Li{x)\ < Cix exp(-C72(logx)2/^-'), 

for any positive e, for certain positive constants Ci and C2 depending on e; this result is due to 
Vinogradov and Korobov in 1958. 

Problem E encodes a modification of a criterion of Guy Robin ||l^ for the Riemann hypothesis. 
Robin's criterion states that the Riemann hypothesis is true if and only if 

cj(n) < e'''nloglogn for all n > 5041, (1-2) 

where 7 ~ 0.57721 is Euler's constant. This criterion is related to the density of primes, as 
explained in §2. Our aim was to obtain a problem statement as elementary as possible, containing 
no undefined constants. However the hard work underlying the equivalence resides in the results 
of Robin stated in §3, where we give a proof of Theorem |1 . 1| . 

Before coming to the proof, in the next section we describe how the Riemann hypothesis is 
related to the sum of divisors function. The connection traces back to Ramanujan's work on highly 
composite numbers, and involves several results of Erdos with coauthors. 



2. Colossally Abundant Numbers 



The Riemann hypothesis is encoded in the criterion of Theorem 1.1 in terms of the very thin set 



of values of o"(n) that are "large." The sum of divisors function is given by 

a{n)= \{{l+p + p' + ...+p'')=n\{{l + - + ... + \), (2.1) 

p°||n p^lln 

where the product is taken over all primes p dividing n and the notation means p"" divides n 

but does not divide n. Most values of (T{n) are on the order of Cn. The average size of o"(n) 
was essentially found by Dirichlet in 1849, and is given by the following result of Bachmann, cf. 
Hardy and Wright 0, Theorem 324]. 
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Theorem 2.1 (Bachmann) The average order of a{n) is \n. More precisely, 



as n oo. 



^(^(j) = T^^^ + 0{nlogn) 



The maximal order of cr{n) is somewhat larger, and was determined by Gronwall in 1913, see 
Hardy and Wright g, Theorem 323, Sect. 18.3 and 22.9]. 

Theorem 2.2 (Gronwall) The asymptotic maximal size of a{n) satisfies 

limsup — - — ^—^ — = e"^, 
n^oo n log log n 

where 7 is Euler's constant. 

This result can be deduced from Mertens' theorem, which asserts that 



e ' 



p log X 

p<x 

as X — > 00. A much more refined version of the asymptotic upper bound, due to Robin Theorem 
2], asserts (unconditionally) that for all n > 3, 

Tl 

cr(n) < e^n log log n + 0.6482- . (2.2) 

log log n 

In §3 we will show that 

7n 



Hn + exp{Hn) log(-ff„) < e^n log log n + 



logn' 

the right side of which is only slightly smaller than that of (|2.2| ). Comparing this bound with ( |2.2D 
shows that the inequality (^]^), if ever false, cannot be false by very much. 

The study of extremal values of functions of the divisors of n is a branch of number theory 
with a long history. Let d{n) count the number of divisors of n (including 1 and n itself). Highly 
composite numbers are those positive integers n such that 

d{n) > d{k) for I <k <n-l. 

Superior highly composite numbers are those positive integers for which there is a positive exponent 
e such that 

d{n) ^ d{k) , „ , , 
——- > for ah k> 1, 

so that they maximize ^^^^ over all n; these form a subset of the highly composite numbers. The 
study of these numbers was initiated by Ramanujan. One can formulate similar extrema for the 
sum of divisors function. Superabundant numbers are those positive integers n such that 

> for 1 < /c < n - 1. 

n k 



3 



Tl 


Factorization of n 


n 




2 








l.tJUU 


6 


2 • 


3 






2.000 


12 


2^ 


• 3 






2.333 


60 


2^ 


• 3 • 


5 




2.800 


120 


2^ 


• 3 • 


5 




3.000 


360 


2^ 


• 32 


• 5 




3.250 


2520 


2^ 


• 32 


•5-7 




3.714 


5040 


2^ 


• 32 


•5-7 




3.838 


55440 


2^ 


• 32 


• 5 • 7 • 11 




4.187 


720720 


2^ 


• 32 


• 5 • 7 • 11 • 13 




4.509 


1441440 


2^ 


• 32 


• 5 • 7 • 11 • 13 




4.581 


4324320 


2^ 


• 33 


•5- 7- 11 -13 




4.699 


21621600 


2^ 


• 33 


• 52 • 7 - 11 • 13 




4.855 


367567200 


2^ 


• 33 


• 52 • 7 - 11 • 13 


• 17 


5.141 


6983776800 


2° 


• 3'' 


• 52 • 7 • 11 • 13 


• 17 • 19 


5.412 


160626866400 


25 


•33 


• 52 • 7 • • • 23 




5.647 


321253732800 


26 


•33 


• 52 • 7 • • • 23 




5.692 


9316358251200 


26 


•33 


• 52 • 7 • • • 29 




5.888 


288807105787200 


26 


•33 


•52 •7- ••31 




6.078 


2021649740510400 


26 


•33 


• 52 • 72 . 11 • • • 


31 


6.187 


6064949221531200 


26 


•34 


.52.72.11... 


31 


6.238 


224403121196654400 


26 


•3^ 


• 52 • 72 • 11 • • • 


37 


6.407 



Table 1: Colossally abundant numbers up to 10 



Colossally abundant numbers are those numbers n for which there is a positive exponent e such 
that 

(Tin) a(k) 

-r-^ > -i-^ for all k>l. 

SO that n attains the maximum value of over all k. The set of colossally abundant numbers 
is infinite. They form a subset of the superabundant numbers, a fact which can be deduced from 
the definition. Table || gives the colossally abundant numbers up to 10^^, as given in ||]. Robin 
showed that if the Riemann hypothesis is false, there will necessarily exist a counterexample to 
the inequality ( |1.2| ) which is a colossally abundant number, cf. []l8| . Proposition 1 of Section 3]; 
the same property can be established for counterexamples to ( |1.1| ). (There could potentially exist 
other counterexamples as well.) 

Superabundant and colossally abundant numbers were studied]^ in detail by Alaoglu and Erdos ^ 
in 1944. As evidenced in the table, colossally abundant numbers consist of a product of all the 
small primes up to some bound, with exponents which are nonincreasing as the prime increases. 
The values of these exponents have a characteristic smooth shape which can be almost completely 
described (see ( |2.3[ ) below). In fact these classes of numbers had been studied earlier, by Ramanu- 
jan, in his 1915 work on highly composite numbers [0]. The notes in Ramanjuan's Collected 
Papers report: "The London Mathematical Society was in some financial difficulty at the time, and 
Ramanujan suppressed part of what he had written in order to save expense." |[l^, p. 339]. Only 
the first 52 of 75 sections were printed. The manuscript of the unpublished part was eventually re- 
discovered among the papers of G. N. Watson, in "Ramanujan's Lost Notebook" (see Andrews |^ 

^Alaoglu and Erdos use a slightly stronger definition of colossally abundant number; they impose the additional 
requirement that ^["l > 'j^ ^'il must hold for 1 < k < n. With their definition the colossally abundant numbers are 
exactly those given by (2.3) below. 
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and Ramanujan [|T5|, pp. 280-308]) and finally published in 1997, in ||T6|. In it superabundant 
and colossally abundant numbers are considered in Section 59, as special cases of the concepts of 
generalised highly composite numbers and superior generalised highly composite numbers, for the 
parameter value s = 1. Ramanujan derived upper and lower bounds for the maximal order of gen- 
eralised highly composite numbers in Section 71, assuming the Riemann hypothesis. (The Riemann 
hypothesis was assumed from Section 40 onward in his paper.) His bounds imply, assuming the 
Riemann hypothesis, that ( |1.2[ ) holds for all sufficiently large n. 

The results of Alaoglu and Erdos in their 1944 paper are unconditional, and mainly concerned 
the exponents of primes occurring in highly composite and superabundant numbers. In consider- 
ing the exponents of primes appearing in colossally abundant numbers, they raised the following 
question, which is still unsolved. 



Conjecture. (Alaoglu and Erdos) If p and q are both primes, is it true that and are both 
rational only if x is an integer? 

This conjecture would follow as a consequence^ of the four exponentials conjecture in transcendental 
number theory, which asserts that if ai,a2 form a pair of complex numbers, linearly independent 
over the rationals Q, and 6i,&2 is another such pair, also linearly independent over the rationals, 
then at least one of the four exponentials {e"'''^ : i, j = 1 or 2} is transcendental, see Lang |l^, pp. 
8-11]. Alaoglu and Erdos ||l]. Theorem 10] showed that for "generic" values of e there is exactly one 
maximizing integer n , and the exponent ap(e) of each prime p in it is 

a,{e) = L^^^M^^ll^Il^M^l^J _ 1; (2.3) 
iogp 

furthermore for all e > this value of n is maximizing. However for a discrete set of e there will be 
more than one maximizing integer. Erdos and Nicolas |^, p. 70] later showed that for a given value 
of e there will be exactly one, two or four integers n that attain the maximum value of ^rr?- 
the conjecture of Alaoglu and Erdos above has a positive answer, then no values of e will have four 
extremal integers. This would imply that the ratio of two consecutive colossally abundant numbers, 
the larger divided by the smaller, will always be a prime, and that every colossally abundant number 
has a factorization ( |2.3| ) for some nonempty open interval of "generic" values of e. 

Now we can explain the relevance of the Riemann hypothesis to the extremal size of a{n). 
Colossally abundant numbers are products of all the small primes raised to powers that are a 
smoothly decreasing function of their size. Fluctuations in the distribution of primes will be reflected 
in fluctuations in the growth rate of taken over the set of colossally abundant numbers. Recall 
that the Riemann hypothesis is equivalent to the assertion that for each e > 0, 

\^{x)-Li{x)\ < 

holds for all sufficiently large x. It is also known that if the Riemann hypothesis is false, then there 
will exist a specific positive constant 6 such that the one-sided inequality 

7r(x) > Li{x) + x^/2+<5 

is true for an infinite set of values with x ^ oo. Heuristically, if we choose a value of x where such an 
excess of primes over Li[x) occurs, and then take a product of all primes up to this bound, choosing 

^For non-rational x consider oi = logp, 02 = logq, 61 = 1 and 62 = a; in the four exponentials conjecture. For 
non-integer rational x a direct argument is used. 
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appropriate exponents, one may hope to construct numbers n with cr(n) exceeding e'^nloglogn by 
a smah amount. If the Riemann hypothesis holds, there is a smaller upper bound for excess of the 
number of primes above Li{x), and from this one can deduce a slighly better upper bound for a{n). 
The analysis of Robin |]l8| gives a quantitative version of this, as formulated in Propositions |3.l| 



and 3.2 below 



One can prove unconditionally that the inequality (1.1) holds for nearly all integers. Even if 



the Riemann hypothesis is false, the set of exceptions to the inequality (1.1) will form a very sparse 



set. Furthermore, if there exists any counterexample to ( |1 . 1[ ) , the value of n will be very large. 
3. Proofs 



The proof of Theorem LI is based on the following two results of Robin |18] 



Proposition 3.1 (Robin) If the Riemann hypothesis is true, then for each n > 5041, 

d < e'^n log log n , (3.1) 

d\n 

where 7 is Euler's constant. 

Proof. This is Theorem 1 of Robin [|^] . ■ 

Proposition 3.2 (Robin) If the Riemann hypothesis is false, then there exist constants 
< /3 < ^ and C > such that 

^y^TX Icj^^ Ic^^^ 

Ed > e"^?! log log n H -5— (3.2) 
(log n)P 

d\n 

holds for infinitely many n. 

Proof. This appears in Proposition 1 of Section 4 of Robin |jl8|. The constant /3 can be chosen 
to take any value 1 — b < (3 < ^, where b = ^{p) for some zero p of C(s) with ^{p) > i, and C > 
must be chosen sufficiently small, depending on p. The proof uses ideas from a result of Nicolas 



1 12 1, Proposition 3], which itself uses a method of Landau. 



We prove two preliminary lemmas. 
Lemma 3.1 For n > 3, 

ex.p{Hn)log{Hn) > e^nloglogn . (3.3) 
Proof. Letting \ t\ denote the integer part of t and {t} the fractional part of t, we have 

^ l<r<t<n ^ r=l 

Thus 

t - {t} {t} 



Hn = l+ I —^dt = logn + 1 - / ^dt. (3.4) 
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From this we obtain 



Hn = logn + j+ I ^dt, 



( 



where we have set 

This is Euler's constant 7 = 0.57721 . . . , since letting n ^ 00 yields 

7 = lim {Hn - logn), 

n^oo 

which is its usual definition. Now ( |3.5| ) gives 

Hn > logn + 7, 

which on exponentiating yields 

exp(i7„) > e^n. ( 
Finally > logn so log{Hn) > log logn > for n > 3. Combining this with ( |3.6D yields (jsT 

Lemma 3.2 For n > 20, 

7n 



Hn + exp(i?„) log{Hn) < e<n log log n + 



logn 



( 



Proof. The formula (3^) implies, for n > 3, 



log Hn < log(log n + 1) < log log n + 



1 



log(n + 1) 



( 



Next, (PD yields 



{t} 

exp{Hn) = exp(logn + 7+ / —^dt 



< e^n exp 



e^n exp 



°° dt 



n 



( 



Since 



we obtain from (3.9) that 



exp(2;) < 1 + 2x for < x < 1 , 



eMHn) < e^n 1 + - 
n 



(3 



Combining this bound with (|3.8| ) yields 

exp(i?„) \og{Hn) < e'n log log n + 
< e^n log log n + 



e'n 



log(n + 1) 
6n 



+ 2e'''(log logn + 1) 



logn 



for n > 10. Using Hn < logn + 1 and > logn + 1 for n > 20 yields (3.7). 



Proof of Theorem 



1.1 



<^ Suppose the Riemann hypothesis is true. Then Proposition 3A and Lemma 3J. together 
give, for n > 5041, 

d < e^n log log n < Hn + exp(//„) log ■ 

d\n 

For 1 < n < 5040 one verifies ( |1.1[ ) directly by computer, the only case of equality being n = 1. 

=^ Suppose (1.1) holds for all n. We argue by contradiction, and suppose that the Riemann 
hypothesis is false. Then Proposition applies. However its lower bound, valid for infinitely 
many n, contradicts the upper bound of Lemma p.2| for sufficiently large n. We conclude that the 
Riemann hypothesis must be true. ■ 
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